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Abstract 



A complete description of the calculation of anomalous dimensions (GLAP split- 
ting functions) is given for parton distributions which appear in space-like pro- 
cesses. The calculation is performed in the light-cone gauge. The results are in 
agreement with the previous results of Furmanski and Petronzio. 



^Permanent address 



1 Introduction 



The light-cone gauge has always occupied at special role in the description of 
hard processes at high energy. It belongs to a class of physical gauges in which 
many of the precepts of the QCD parton model are true, because in this gauge 
collinear divergences occur in diagrams corresponding to the parton cascade. We 
thus retain the probabilistic interpretation of a hard scattering event, which is 
obscured in covariant gauges. In fact, by introducing an additional gauge vector 
we obtain many of the advantages of an infinite momentum frame formulation 
in a covariant notation. 

However, doubts have been raised about the utility of light-cone gauge in 
practical calculations One refutation of these misgivings is provided by the 
classic calculation of the two loop splitting functions or anomalous dimensions 
as given by Curci, Furmanski and PetronzioQ for the non-singlet case and the 
calculation of Furmanski and Petronzio[Q for the singlet case. The calculation of 
ref. [Q, however, has never been fully documented. It could be be that this lack 
of complete documentation has acted as a barrier to further developments along 
this line. One example of a further application of this method is the calculation 
of the polarized two loop splitting functions, recently presented in ref. 0. It 
therefore seemed a valuable addition to the literature to provide a more complete 
description of the calculation of the singlet evolution probabilities. This is the 
modest aim of this paper. 

In our calculation n is a light-like vector which serves to define the longitudi- 
nal direction. The momentum of the incoming parton (taken to be massless) is 
denoted by p. Thus we have two light-like vectors which are defined such that, 

71^ = = 0, n ■ p = pn 7^0, n-i = p- < = 0, (1) 

where t is any vector in the transverse plane. In addition to specifying the 
longitudinal direction we also use the vector n to fix the light-cone gauge: 

71 • A = . (2) 

Following reference |^ we use the principal value (PV) prescription to regulate 
the divergences which occur in the light-cone gauge propagator in loop and 
phase space integrals, i.e. 

(3) 




{n ■ kY + S'^{pn) 



This prescription is at variance with the Mandelstam-Leibrandt (ML) treat- 
ment of the l/{n-k) singularities . The ML prescription, since it permits the 
Wick rotation in virtual diagrams, leads to power counting rules and a proof of 
renormalizability of the theory in this gauge. Indeed the one loop non-singlet 
splitting function has been investigated in this gaugelq]. We have chosen not 
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to follow the ML prescription. The technical reason is that it leads to a prolif- 
eration of graphs, because of ghost-like contributions associated with the n ■ k 
propagator^. The physical reason is that we wanted to stay as close as possi- 
ble to old-fashioned perturbation theory in which manifest Lorentz covariance is 
sacrificed in order to have a simple form for unitarity. The advantage of Eq. (||) 
is that the unitarity of the theory is explicit. In this way we hoped to gain a 
greater physical understanding of the two loop anomalous dimensions. It would 
be interesting to repeat the calculation of the two loop anomalous dimensions 
using the ML prescription. 

Of course, the calculation of any gauge invariant quantity such as the two 
loop splitting function is independent of the gauge in which the calculation is 
performed and of the method of calculation. However, the discovery of compli- 
cations in the covariant gauge calculations!^,^ makes the method outlined in 
ref. [0 and implemented in refs. even more attractive. Not only is this 
method close to the parton model, but it also leads to compact answers and 
may be the most efficient method from a calculational point of view. It might 
present a viable method for the analytic calculation of the three loop splitting 
function^. 



2 Calculation of anomalous dimensions 



2.1 Factorization 

In this section we shall explain the method of factorization of the two-particle 
irreducible (2PI) diagram in the light-cone gauge. It is not our intention to 
repeat the discussion which is clearly provided in ref. |^ . We only include those 
details which are necessary to present the structure of the calculation or to 
define the notation. Following ref. we define a generalized ladder expansion 
by introducing the 2PI kernel Kq: 



M = Coil + Ko + K^+K-^ + ...) 



Co 



l-Ko 

Factorization occurs by introducing the projector onto physical states, V, 
1 1 



(4) 



1 



Ko 1 - (1 - V)Kq - VKo 

_ ' 1 

l-{l-V)Ko 

Defining the modified kernel K, 

K 



\-VKo[l-{l-V)Ko\ 



Ko 



1 - (1 - V)Ko 



(5) 



(6) 



^Results for low moments of the three loop splitting functions are given in Ref. 
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we can thus write M as 



= c r (7) 

where 

^ = ^°l-(l-mo ' 

r = — i — . (9) 

At this stage the factorized structure becomes apparent. In the hght-cone gauge 
the 2PI kernels Kq are finite before the integration over the sides of the ladder 
is performed. CoUinear singularities appear only after integrating over the lines 
connecting the rungs of the ladders |p^. All coUinear singularities are contained 
in r, whereas C is interpreted as the (finite) short distance cross section. Re- 
expanding we find that 

= i^o (1 + (1 - V)K^ + (1 - V){K^{1 - V)Ko) + ...), (10) 



r = 1 + VK+{VK){VK) + ... 



(11) 



Restoring the indices and regulating coUinear singularities by going to d = 4 — 2e 
dimensions we have that 



1 + PKo + V{Kl) - P{KoPKo) 



and explicitly 



Tij{z, Us, -) Zj 



df^k n -k 

S(z )UiK 



(27r 



pn 



l-PK ^ 



(12) 



(13) 



where 'PP' extracts the pole part of the expression on its right and Zj {j = q{g)) 
is the residue of the pole of the full quark (gluon) propagator, contributing to the 
diagonal splitting functions. Furthermore the spin averaged projection operators 
onto physical states are given by, 



d~2 



pn 



(14) 
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2.2 Derivation of GLAP equation 

The property of factorization allows us to separate the low momentum physics 
from the high momentum physics in a multiplicative way. This separation is 
performed at a scale ^, which is completely arbitrary, and no physical prediction 
can depend on it. In this section we investigate the constraints provided by this 
condition. For simplicity, we will consider a non-singlet cross section which can 
only be initiated by a quark. We therefore have the factorized result. 



— , e) = — , as(M^)) ® r,g(as(M^), e) (15) 



where we have indicated that Tgq does not depend on (i.e. Q^/iJ.^), which is 
a consequence^ of the finiteness of the kernel Kq in the light-cone gauge. The 
symbol (8) indicates a convolution integral over longitudinal momentum fractions 
of the type 

f®g=l dy dz f{y)g{z) 5{x -yz) . (16) 
Jo 

If we take moments, 

r-l 

7-1 



fU) = / dxx^-'f{x) (17) 
on both sides of Eq. (p^), it reduces to a simple product: 

^(j, ^,as(M^),e) = ^q{j, ^,a5(/^^))rgg(j,as(M^),e) • (18) 

ffq is the short distance cross section from which all singularities have been 
factorized. Tgq contains the mass singularities which manifest themselves as 
poles in e. The independence of the full cross section of /i implies that 

:a = (19) 



din ^■ 



and hence that 

—^\nTqq{j,as{^l'^),e) = -—^\naq{j,^,as{^i'^)) = lqq{j,as{^i'^)) 



(20) 

The function is known as the anomalous dimension, because it measures the 
deviation of dq from its naive scaling dimension. It must be finite and can only 
depend on as(/x^) because these are the only variables common to both Tqq and 
dq. The anomalous dimension is extracted from Eq. ( pO| ) in the following way: 
Because the /i dependence of Tqq enters only through the running coupling we 
have that, 

Iqqij^asiii^)) =f3{as,e)- — \nTqq{j,as{fJ.'^),e) , (21) 

das 
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where (3{as, e) is the d-dimensional QCD (3 function in the MS scheme, 



P{as,e) = -^^ = -eas + P{as) . (22) 



In the MS scheme Tqq is given by a series of the form 



r9gU,as,e) = 1 + 2^ ^ ■ (23J 



Comparing the coefficient of the term of order e° we find that 



Integrating Eqs. (|20|,pi|) one obtains 



{rots 
J 



and 



asifJ-"^), e) = 1, a5((3^)) exp < J dX 



/3(A) - eA 



(26) 



with the running couphng as{Q^)- In order to obtain the hadronic cross section, 
a{Q'^/fi'^, as,e) has to be convoluted with 'bare' ('unrenormaHzed') quark den- 
sities q{as,£) which contain mass singularities that must exactly cancel those 
in Tqq. The resulting 'dressed' ('renormalized') quark distribution function 



qU,Q^) = exp<^ dX 



>qU,as{p'),e) (27) 



/3(A) - eA 

is free of mass singularities and satisfies the non-singict evolution equation 
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lggU,as{Q'))q{3,Q') (28) 



dlnQ2 
which, defining 

/ dzz^~'^Pqq{z,as)=lqq{i,as), (29) 
Jo 

in X space takes the form of the GLAP equation [p^ 13 1: 

^1^^ - dyj'^ dz Pqq{y, as{Q'))q{z, Q^) S{x - yz) . (30) 
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as) = {^) P(°)(.) + (If + . . . (31) 



(32) 



Expanding 
one has 

T,,{x, as, e) = <5(1 - x) - 1 (^^P^^K^) + i [^f pW^^) + . . 

The generahzation to the singlet case is straightforward. 
2.3 Non-singlet and singlet equations 

The separation into singlet and non-singlet parts depends on the properties of 
the kernel. Using SU (/) flavour symmetry we may define the following combi- 
nations of qq and qq matrix elements: 

p s pV I pS 

P - — t) , P^ 4- P''' 

P± = P^±P^. (33) 
In addition, because of charge conjugation invariance, we have that 



p 


= P- - 




P - 


— P- 




P 


— P- 




p 


= p - 
^gqi ■ 


(34) 



At two loop order, there is a non-zero contribution from and P^, but 
we have the additional relation 

Pqq = Pqq ■ (35) 

which simplifies the treatment of the non-singlet pieces. 

For each flavour we define the sum and difference of the quark and anti-quark 
distributions as 

qf =qi±qi. (36) 
One then finds that the combinations 

V, = qi (37) 

and 

fe 

Ti=Y,<lt-kqt (38) 

i=l 
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(where i,k — 1, . . . / = fc^ — 1) arc non-singlets, i.e., evolve according to 

Eq. ( pO[ ) with the kernels and P+, respectively. 
The singlet Altarelli-Parisi equation ifQ[^2[[l3| 



d\nQ^\ G{j,Q^) 



S(j,Q2) 



(39) 



where G(j) is the moment of the giuon distribution and is the singlet quark 
combination. 



/ / 



(40) 



The elements of the anomalous dimension matrix are given in terms of the 
kernels defined in Eqs. ( ^3[]35| ) as, 



Pa. 



P^ + nfiP, 

'^nfPq^g 



qq' 



P — P 
^gq ~ ^gqi 



(41) 



2.4 Renormalization constants 

The notation for the renormalization constants is shown in Fig. |^. We define 
the integral 

which contains the divergences in the PV regulator S (see Eq. (|)) arising from 
the light-cone gauge propagator. As already noted in ref. [D, use of the PV 
prescription (||) in the light-cone gauge entails the disagreeable feature that the 
renormalization constants depend on the longitudinal momentum fractions x. 



Zg{x) 
Zg{x) 

Z'^^\xi,X2,X^) 
Z^^\xi,X2,Xz) 



Uq 1 

1 + — — 

2tt 2e 



= 1 



1 + 



as 1 

as 1 
2tt 2e 



1 + ^1 
2tt 2e 



CF(-3 + 4/o + 41na;) 



ATt 11 

^ + A^c(-y +4/o + 41na;) 



Cf(3- 4/o - 21n2;i - 21na;2) - 2iVc(/o + Inxg) 



11 ATt 
Nc{^ - 6/o - 2 Inxi - 2 ln.X2 - 2 In 0:3) - ^ 



(43) 



'Note that the notation for the off-diagonal terms is different than in ref. 
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Zg(x) 



7 '^V ^ 



X ' Xl 

Figure 1: Renormalization constants and the vertices which they renormahze 



where 



Nc = 3, T/ = TRUf 



(44) 



When combined in the appropriate combinations to investigate the renormahza- 
tion of the charge, the momentum dependent terms and the divergent integrals 
/o canceL Thus the relationship between the bare and renormalized couplings 
is 



-,(0) 



aS^^^Z^^^ {Xl,X2,X3)^ Zg{Xl)Zg{X2)Zg{X3) 

as 1 , lliVc 2Tf 
27r 2e ^ 6 3 



(45) 



asfJ- 



as^l 



as 1 
27r 4e 



/3o 



2.5 Topologies of NLO graphs 

The basic topologies of all 2PI diagrams which occur in two loops are shown 
in Fig. 1^. The notation of the topologies (b)-(i) is determined by the labelling 
of the diagrams for the non-singlet calculation given in ref. We have not 
included diagrams which can be obtained by reflection about the vertical axis 
which occur in cases (c),(d),(e) and (j). Topologies (hi) correspond to the terms 
ViK'^) - ViKoPKo) in Eq. (|l2|), aU other topologies belong to VKq. As an 
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example, the diagrams corresponding to are given explicitly in Fig. ^. Fig. ^ 
shows the diagrams for (b) and for P^ (h,i). The appendices give the 
necessary ingredients needed for the evaluation of the real and virtual graphs in 
Fig. |. 




(b) (cd) (e) (fg) 




Figure 2: Basic topologies of the diagrams 



3 Results 

As in Eq. ( |3l| ) we define the perturbative expansion 

i^.(^,..)=(g)/f)(.)+(|f)^if)(.) + ... (46) 
The full one loop results are included for completeness [|2| 

x)} (47) 
(48) 
(49) 

x)} + ^<5(l-x). (50) 



P(°)(x) ^ 2T,{x' + {l-xf} 
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(b) (c) (d) (e) 




Figure 3: Diagrams for the qq part of the NLO non-singlet splitting functions 

In order to write the full result for the two loop splitting functions P^^^ we 
introduce the notation 

2 

Pqq{x) = -l-X (51) 

Pqg{x) = x2 + (1-.t)2 (52) 
Pmix) = ^^^^^ (53) 

Pgg{x) = —^ + --2 + x{l-x). (54) 
1 ~ X X 
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CfTr 


Terms 


(hi) 


X ln^{x) 


-1 


\n^{x) 


-1 


x"^ ln(x) 


8/3 


X ln(x) 


5 


ln(x) 


1 




-56/9 


X 


6 


1 


-2 


l/x 


20/9 



Table 1: Results for the gg-Singlet diagram 



The graph-by-graph results are given in Tables [T]-4, where we only list the 
contributions to the singlet splitting functions since those for the non-singlet 
case have been presented in a similar table in ref. Q. Our final full results 
for the two loop non-singlet splitting functions read for a; ^ 1 : 



^9?'^ = Cl[ - [21nxln(l - x) + - In x] p,, (a:) 

3 7 1 1 

^(- + -x) Inx — -(1 -|- a;) In a: — 5(1 — a;)| 

+CfNc[[]^ \n x+^\nx + ^- y ]pg,(x) + (1 + x) In a; + y (1 - x) } 
+CpTf[ - Inx + y]p,,(a;) - ^(1 - a;))} (55) 
Pj^'^'^ =CF(CF-^){2p,,(-a:)52(a;) + 2(l + a;)lna; + 4(l-a;)} . (56) 
Our results for the singlet terms areQ, 

(57) 

+2CFTf{— - 2 + &X- —x^ + (1 -h5a;+ -a;2)lna;- (1 -Fa;)ln^T| 
L 9a; 9 3 J 

^ CpTf^A - 9x - (1 - Ax) In a; - (1 - 2x) \t? x + 41n(l - x) 
+ [21n' (^) - 41n(^) - + 10]p,,(a:)} 
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_ fl82 14 40 ,136 38,, , , ^ „ x , 2 

+^cr/|— + -9-2; + — + i—x - — ) Ina; - 41n(l - .x) - (2 + 8x) In^ x 

+ [ - \n^x+^\nx - 21n^(l - a:) +41n(l - x) + ^ - '^]pqg{x) 
+2pgg{-x)S2ix)} (58) 

Pgl^ = C|{ - ^ - ^a; + (2 + Ix) Inx - (1 - ^a;) In^ x - 2a:ln(l - x) 
-[31n(l - x) + ln'(l - x)]pgg{x)} 

+CfNc\^ + ^x + ^x^ - (12 + 5a;+ |.T2)lna;+ (4 + x) ln^T + 2xln(l 
L 9 lo 9 o 

+ [ - 21nxln(l -x) + i In^ a + y ln(l - a) + 111^(1 - x) - y + ^jPg-jW 
+52(x)pgg(-x)| 

+CFTf[ - - [y + ^ ln(l - x)]pg,(a^)} 



p(i) = CfT/I - 16 + 8x + + 7^ - (6 + lOa) In a; - 2(1 + a;) In^ x] 
L 3 3a J 

+iVcr/{2 - 2a + y (a:2 _ i) - 1(1 + .x) In a - yP<,<,(a:)} 

67 TT^ 1 

+ [— - 41na;ln(l - x) + In^a - y]pg3(a) + 2pgg(-.T)S'2(a))| (60) 
where the function 6*2 (a) is defined 

52(x)= _ln(^). (61) 

In the small-x limit S2 becomes 

^2-lln'x-y +0(x) . (62) 

All results in Eqs. ( p5[]60| ) are in complete agreement with the corresponding 
results in . They can be extended to all values of a using a trick to evaluate 
the endpoint contributions in Eqs. ( |55| , |60| ). The sum rule from the conservation 
of fermion number is 

dx (^Pqq{x) - Pgg{x)^ = j dxP-{x)^0. (63) 



*Note that the definition of 52 in ref. M contains a typographical mistake. 
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The conservation of momentum leads to the foUowing two relations, 



dx x(^Pgq{x) + Pgqix)j = , 

J dx x(^Pgg{x) + Pgg{x)^ = . 



(64) 
(65) 



These results for the integrals of the splitting functions are satisfied if one makes 
the substitutions 

in Eqs. ( pl| , |5^ ) and adds in the end-point contributions to Eqs. ( |55| , |60| ) |Q, 



p^.(l) ^ pV'W _|_ 
99 99 ^ 



c 



17 Utt^ 



18 



(5(1 -a;) 



p(i) ^ p(i) 



A^c{3+3C(3)}-Cf^T/--7VcT/ 



3C(3)} 
(67) 

6(1 - x) (68) 



where C(3) ~ 1.202057. The substitution in Eq. ( |66| ) is obviously not necessary 
if the factor of 1/(1 — a;) has a coefficient which vanishes a.t x — 1, such as ln(x). 
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CpTf 




Terms 


(cd) 


(e) 


(fg) 


(hi) 


Sum 


(b) 


(cd) 


(hi) 


Sum 


Pqg{x) ln^(l — x) 


2 




-2 


2 


2 






-2 


-2 


Pqg{x) ln'^{x) 


2 


-4 




4 


2 


-1 






-1 


Pqg{x) ln(a;) ln(l — a;) 


4 


-8 


-4 


4 


-4 




4 


-4 




Pqg(x) lo (ln(l - x) + \n(x)) 


8 


-8 


-4 


4 






4 


-4 




Pqg{x) \tl{x) 


-4 


8 


3 


-3 


4 


3 




35/3 


44/3 


Pqg{x) ln(l — x) 


-2 


6 




-8 


-4 


4 






4 


Pqg{x) Li2(l - X) 


-8 






8 






8 


-8 




Pqg{x) TT^/S 


6 


-4 


-2 


-2 


-2 




-1 


2 


1 


Pqg{-X) S'J 












2 






2 


Pqg(x) h 




8 




-8 




4 


-4 






Pqg{x) h 


-8 


8 


4 


-4 






-4 


4 




Pqg \X ) 


-14 


8 


7 


9 


10 


-13 


-4 


-65/9 


-218/9 


X \ii^(x) 








2 


2 






-8 


-8 


lii^ix) 








-1 


-1 






-2 


-2 


X lii{x) 








4 


4 


20 


10 


46/3 


136/3 


lii{x) 


2 


-8 




5 


-1 


1 


-4 


-29/3 


-38/3 


ln(l - x) 


4 




-4 


4 


4 






-4 


-4 


lo 


8 


-8 


-4 


4 






4 


-4 




X 


2 




1 


-12 


-9 


-4 


4 


14/9 


14/9 


1 


-8 


6 


3 


3 


4 


17 


2 


11/9 


182/9 


1/x 
















40/9 


40/9 



Table 2: Results for the qg diagrams 









CpTf 


Terms 


(cd) 


(fg) 


(hi) 


Sum 


(b) 


(cd) 


(e) 


(hi) 


Sum 


(e) 


Pgq(x) ln^(l - x) 


1 


-1 


-1 


-1 








1 


1 




Pgq{x) lll^{x) 










-1/2 


1 


-2 


2 


1/2 




Pgq{x) ln(a;) ln(l — a;) 


4 


-2 


-2 








-4 


2 


-2 




Pqq(x) lo (Infl - x) + \n(x)) 


4 


-2 


-2 






2 


-4 


2 






Pgq(x) \n(x) 


-3 


3/2 


3/2 




-3/2 






3/2 






Pnn(x) ln(l — X) 

i gq \ / V / 


1 


-2 


-2 


-3 


2 




11/3 


-2 


11/3 


-4/3 


Pnn(x) Li2(l — X) 


4 




-4 






-4 




4 






Pnn(x) 7r^/3 

I / 1 




-1 


1 






5/2 


-2 


-1 


-1/2 




Pgq{-x) S2 










1 








1 




Pgq{x) lo 


4 


-2 


-2 




2 






-2 






Pgq{x) h 


-4 


2 


2 






-2 


4 


-2 






Pgq{x) 


-10 


5 


5 




-13/2 


-2 


67/9 


14/9 


1/2 


-20/9 


X \ii^(x) 






1/2 


1/2 








1 


1 




lii^ix) 






-1 


-1 








4 


4 




x^ ln{x) 
















-8/3 


-8/3 




X lii{x) 


1 




5/2 


7/2 


-1/2 


4 


-4 


-9/2 


-5 




lii{x) 






2 


2 


-10 


-5 




3 


-12 




X ln(l — x) 


2 


-2 


-2 


-2 








2 


2 




X lo 


4 


-2 


-2 






2 


-4 


2 






x' 
















44/9 


44/9 




X 


-4 


3/2 


-1 


-7/2 


17/2 


1 


11/3 


-86/9 


65/18 


-4/3 


1 


1 


1/2 


-4 


-5/2 


-2 


2 




28/9 


28/9 





Table 3: Results for the gq diagrams 





CpTf 


NaT, 




Terms 


(b) 


(hi) 


Sum 


(b) 


(rd) 




(fg) 


Sum 


(b) 


(cd) 


(e) 


(fg) 


(hi) 


(jk) 


Sum 


Pgg{x) \n{l - x) 




















2 




-2 








Pg^(x) \n (x) 


















-1 


2 


-4 




4 




1 


Pgg(x) ln(.T:)ln(l - x) 




















8 


-8 


-4 






-4 






















12 


-8 


-4 








Pgg{x) \n{x) 










4/3 




-4/3 






-11/3 




11/3 








Pgg(x) ln(1 - x) 










8/3 


-8/3 






4 


2/3 


22/3 


-4 


-8 






Pgg(x) 7rV3 























A 

-4 


-I 






1 
- 1 


Pgsi *J '-'2 


















2 












2 


Pgs 


















4 


8 




-4 


-8 






Pgg(^) h 




















-12 


8 


4 








Pgg(^) 










40/9 


-40/9 


-20/9 


-20/9 




-170/9 


134/9 


103/9 






67/9 






-2 


-2 




















4 




4 






-2 


-2 




















4 




4 


\n{x) 


8/3 


-8/3 




-4/3 


4/3 








-11/3 


-11/3 






-22/3 




-44/3 


x \n{x) 


8 


-18 


-10 


-4 


8/3 






-4/3 


-17/2 


-23/6 
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11/3 


1ti(.t) 


8 


-14 


-6 


-4 


8/3 






-4/3 


-17/2 


-23/6 






4 




-25/3 


ln(2;) /x 


8/3 


-8/3 




-4/3 


4/3 








-11/3 


-11/3 






22/3 






x-^ 


-76/9 


136/9 


20/3 


38/9 


-8/3 




4/3 


26/9 


136/9 


22/3 




1/3 


-46/3 




67/9 


X 


4 


4 


8 


-2 


5/3 




-5/3 


-2 


-105/8 


-25/3 




-5/12 


19/2 


-9/8 


-27/2 


1 


-4 


-12 


-16 


2 


-5/3 




5/3 


2 


105/8 


25/3 




5/12 


-19/2 


9/8 


27/2 


\/x 


76/9 


-64/9 


4/3 


-38/9 


8/3 




-4/3 


-26/9 


-136/9 


-22/3 




-1/3 


46/3 




-67/9 



Table 4: Results for 



the gg diagrams 



4 Summary 



This paper has presented a recalculation of the two loop anomalous dimension 
for space-like processes. The results given in this paper are therefore not new. 
The new features arc the presentation of the results in a coherent notation, the 
description of some of the integrals which are required to derive the results and 
a detailed description of the contributions of the sub-diagrams to the results. 



A Virtual Integrals 

A.l Two point function 

The evaluation of the virtual integrals involving non-covariant denominators of 
the form l/{n-k) requires some care. We define 

Z+ = n-l 
I- = p-l 

dl^l = dl+dl-d'^-Hi. . (A.l) 

We shall evaluate the integrals by explicitly performing the integrals over l±_ 
keeping 1+ fixed. This formulation will be useful when /(/+) contains poles in 
1/1+ coming from non-covariant denominators as long as the method used to 
regulate the 1+ singularity does not involve I- . The general two-point function 
then reads 



/ 



d'^l /(/^ 



{2-kY {P + ie){{l - kf + ie) 
At: y r(l + e) 



167r2 V -k^ 



e 



dz fil+) z-'il-z)-' (A.2) 



where / is an arbitrary function, d = 4 — 2e and z = l+/k+ the boost invariant 
rescaled value of Z+. If /(/+) = 1 we recover the normal covariant result, 

^ _ , d'^l 1 i / 477 V r(l + e)r2(l-e)l 



r dH 

J W) 



tY {P + is){{l - kf + ie) \-k^ J r(2 - 2e) e' 

(A.3) 

The result for the integral with one non-covariant denominator is 



Y \l+ - k+ j {P + ie){{l - kf + ie) 
47r Vr(l-he)r2(l-e) f /o + lna; 



167r2 V-fcV r(l-2e) | e 
-Io\nx + h-Un^x-^+0{e)\ (A.4) 
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where x — n ■ k/pn and we have indicated that the PV prescription defined by 
Eq. (^) has been used. Jq is defined as in Eq. (^); in the small-i5 Umit it reduces 
to 



Furthermore, 



For the case without an endpoint singularity in the integral over the plus 
component we may take the limit (5 — > and hence obtain, 

k+ 1 



{2TT)'i 1+ ~ p+ {P + ie){{l - A:)2 + ie) 

47r y r(l + e)r2(l -e) f ln(l - x) 



16.H-^V r(l-26) I e ■2Li,(l-x)-- 
+21nx-ln(l -x)-^ ln^(l - x) + 0(e)| (A.7) 

where Li2(x) is the usual dilogarithm fimction, 

U2{x)^- H-^^^^—^dt. (A.8) 
Jo t 

We note that one also needs the vector two point function with one non-covariant 
denominator, 

Assuming Lorentz covariance of the momentum integral one finds 

= [k^ - ^n^) J. + ^n^' {-J... + ^ Jr,„) (A.10) 

where 

d'^l „,/ k+\ 1 



It turns out that the integral Ji^„ always cancels out in the final answer. 
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(a) 



P3 



p-k 




Figure 5: (a) Vertex correction graph (b) One parton emission (c) Two parton 
emission. 



A. 2 Three point function 

We shall only consider the special case which is needed for our purpose. We 
employ the momentum assignments ^ {p — k)^ = and define the boost in- 
variant quantities, x = k+/p^, y — l+jp^, z ~ yjx ~ l^/k+. The corresponding 
diagram is shown in Fig. One finds 



-i / 47r y r(l + e) 



{2TTy {P + ie){{l - k)^ + ie){{l - p)^ + ie) Wn^k^ \-k^ 
dyf{l+)z-\l-z)-'-^ - ,.<l-x) 



z-l 



(A.12) 



where 2F1 is the hypergeometric function. With a little work one can show that 
for the special case f{l+) = 1 one recovers the normal covariant result, 



{2ttY [P + ie){{l - ky + ie){{l - pf + le) 
i ( 47r \'r(l + e)r2(l-e) 



167r2fc2 



r(l-2e) 



(A.13) 



The explicit result for the three point function with one light-cone gauge denom- 
inator using the PV prescription is (we have performed a shift of the integration 
variables relative to Eq. (A.12)): 



dH 
(2^ 



PV 



1 



1++P+J {{1+ P? + i^W + + P'Y + is) 
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47r y r(l + e)r2(l - e) j 1 Inx-/, 







IGir'^k^ \-ky r(l-2e) 



^2 



+/i-/olna;-2Li2(l-a;)-iln'x-y + 0(e)| . (A.14) 

The other integral which we need can be obtained from Eq. ( |A.14 ) by exchange 
of p and p' 



{2TT)d- \l++p'^J ((;+p)2+ie)(p + ie)((;+p')2+j£) 
47r y r(l + e)r2(l -e) f 1 Ina; - 21n(l - x) - /q 



167r2fc2 V-A:^ y r(l - 2e) 



e 



1 Stt^ 

+/i - /olna; + 2Li2(l - x) - - hi^ a; + In^Cl - x) + 0(e) ^ .(A.15) 

2 6 I 

A useful relation in comparing these results with the real diagrams is 

r(l + e)r(l-e) = l + e2!L! + 0(e4) _ (a.i6) 

6 

B Real integrals 

In this appendix we will describe some of the integrals which occur in diagrams 
involving the emission of real partons. As illustrated in Figs. ^(b,c), we will 
denote the momenta of the emitted particles by p^ and p4 and the momentum 
of the 'observed' parton line by k. The phase space for one and two parton 
emission, and keeping k'^ and n ■ k fixed, is given by (we set pn = 1 in the 
following) 



PS 



(1) 



d'^k 5{x~n- k) 5{\k^\ + {p- P7if)5'^{p - Ps - k) (B.l) 



J d^k 6{x~n- k) d{\k^\ + {p-p3- Pif)5''{p -P3-P4-k). 

(B.2) 
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Integrating over irrelevant angles in d dimensions we have for the transverse 
phase space, 



If the integrand is independent of 6i we can integrate further to obtain 



d'^'^kr - — r / dk^ k^"^' . (B.4) 



B.l Crossed ladder diagrams: topology b 

Here we shall describe the integrals needed for the evaluation of the crossed 
ladder diagram (topology b) as shown in Fig. || We introduce a notation for 
the real parton momenta such that, 

Zzi 

Pt = Z2P^ + ^n^" -t2'' (B.5) 

ZZ2 

with transverse momenta ti, t2. In terms of these variables the denominators 
which occur in the diagram in Fig. ^ can be written as 

pI = , pI = , k^ = [p-ps. -Pif = -aiti^ -a2t2^ - 2ti • t2 (B.6) 

where 

(1 - Z2) (l-2l) 

ai = , 02 = . (B.7) 

Zl Z2 

The general form of the matrix element which has to be integrated over the 
phase space of Eq. (|B.2|) is 



A(zi, Z2) + B{z^,Z2)^-^ + C{z^,Z2)^-^ + D{z,, Z2) ^%lf . (B.8) 

The integrations over ti and t2 are finite at small transverse momenta, so that 
before the fc^ integration the expression is finite. Introducing the constants 

F = ^4")^ 



167r2r(l -e) 
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Figure 6: An example of the graph of topology (b) 



we have in the frame specified by Eq. ( ]B.5| ): 



Zi Z2 if 



(B.IO) 



Since the integrals over the transverse momenta are finite we may take the limit 
d — > 4. The values of the integrals for th e inte grands which occur in Eq. (B 



and the phase space weight given in Eq. ( B.f ) are collected in Table BT. The 
remaining one dimensional integrals are easily performed. 



B.2 Real diagrams: topology cd 

The calculation of the topology (b) graphs was facilitated by the fact that before 
integration over the kernel was finite. The situation is more complicated for 
topology (cd), which has two cuts involving either one or two real partons. The 
cut with two real partons (gluons) is shown in Fig. ^ which also shows the 
definition of the kinematics. 



22 



Integrand 


Integral in units of F"^^-^ J dzidz2 S{1 — zi — Z2 — x)9{zi)9{z2) 


< 1 > 


1 


< ^ > 


21 
1-22 


^ ti.t2 \ 

^ 


22 
1-21 


. (tl.t2)" . 


I 2 ZiZ2 (l-2i)(l-22) y J 



Table 5: Real integrals for crossed ladders 



We perform a light-cone decomposition of the light-like momenta p3 and p^: 

pt^ = (i_z)pM + _ii_^M_ti'' 

p'^ = z(l _ y^A' + _JL_„A' _ (B.ll) 

where zy = x. It is expedient to perform a change of variables, 

-ri 



1 


— z 




- y 


1 


- y 




— z 



t2 ^ \/^(r2-(l-z)ri) , (B.12) 
so that p3 and p^ become 



With this choice the propagators of the diagram in Fig. ^ can be written as 

\pI\ =-{p-Pz? =T^ (B.15) 

1 -y 
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Figure 7: An example of the graph of topology (cd) 



\pI\ 



n ■ p3 = 1 - z, 
Hence we have that in this frame, 



J {l~z) (1-2/) ^ ^ ' 

s(\k'\ = 



rf + r^- 2ri • r2 

z 

n ■ p4 = z(l - y) . 



(B.16) 
(B.17) 



/ / 



(B.18) 



If we are integrating over quantities which do not depend on angles we may 
write Eq. (B.18) in the form 



p5(2) ^^2 1^2,1-2. I'dz 



X X \{'^- Z){Z-X) 



^ ^ du uj-\l - uj)-" (B.19) 



where the rescaled transverse momentum lo is defined as 



\P\\ 



{l-y) X 
The results for the integrals are given in Table |. 



(B.20) 
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Integrand 


Value of integral in units of fc^ ^'^ pj^^^j 


< 1 > 


Ifc^l rl dz 

1 — 2€ Jx X 




e 


|p?l 


1 fl dz 1 

€ J X Z \ 


^(l-z)(2:-x)] 


e 




1 fl dz 
e Jx (1- 


x) ^{l-z)(z- 


^)) 


^ \pi\\pi\ > 


2 1 f 1 


{l-z)W. - 





Table 6: Real integrals 



We note at this point that Eq. (B.19) is also suitable for dealing with a 
light-cone gauge denominator term like l/{n ■ p^) in the matrix element. For 
instance, we obtain 



1 



^2|^2|i-2.1 f\^pY(j—] ^F^\P\'-^'-iIo + Hl-x 



< >^ F^lk^l'-^'- I dzPVI ^ I =i^^|fc^|""^'-(/n + ln(l-x)) 

(B.21) 

where Iq is as defined in Eq. ([42|). 

If we have a denominator which depends on the angle between ri and r2, 
the integral is more complicated. For example when we have the denominator 



as given by Eq. (B.16) the angular integration splits into two regions, > 
and rl > rf: 



< 



dz 



(1 — z) X \ (1 — z){z — x) 



I{a) 



(B.22) 



/(a) 



L 


l + a 

( 












4 



auj ^ 



a J ji— OLUj \ 

where 

a = ^ . (B.24) 

x(\ — z) 
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Now by redefinition of variables the two integrals give 



where 



I{a) = J{a) + -J{-) 
a a 



J{a) = / d?;t;-'(l + ai;)"^+^'2i^i(l,l + e;l-e;w) 

^ -^(y-^l + a)-i+2Si^i(l, -2e; 1 - e; 

/£ 1 + a 



(B.25) 



(B.26) 



Hence combining using the identity, 
2Fi(l, -2e; 1 - e; z) + 2^1(1, -2e; 1 - e; 1 - z) = 2 



r^(i-e) 

r(l - 2e) 



one obtains 



/(a) 



ir2(l-e) 1 



z\l~zy (B.27) 



(B.28) 



er(l-2e)l + a 

Thus the final result is as given in Table |. This result can be obtained much 
more easily by performing a shift of the transverse momenta so that only 
depends on a single transverse momentum. 

If we now add a second denominator such a shift is no longer useful. The 
scalar integral with two denominators is given by a simple modification of 
Eq. (IIH), 



< 



\pIM 



> 



p2 |fc2|-l-2e 



dz 



l+2e 



a 



1 — Ix) 



duj - t^)-^2i^i(l, 1 + e; 1 - e; ^-^) 



(B.29) 



where a is as given in Eq. (B.24). By change of variables this integral may be 
further written as 



< 



1 



b?l|pil 



where 



K{a) 



> = F2|fc2|-l-2e 



dv 



dz 



■K{a) (B.30) 



1^-^^+^) + 1) 2^^i(i,i + e;i-e;^^) 

^ dww-^ (^+w)^%Fi(l,l + e;l-e;t;)| . (B.31) 
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The partial result for K(a) which is sufficient for our purposes is 



K{a) 



2r^(l-e) 
er(l-2e) 



( 



1 + eln 



( 



z(l — x) 



Z — X 



) 




So the final result for the integral is as given in Table |6[ 
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